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ABSTRACT 

Our goals are (i) to search for BAO and large-scale structure in current QSO survey 
data and (ii) to use these and simulation/forecast results to assess the science case 
for a new, > 10 x larger, QSO survey. We first combine the SDSS, 2QZ and 2SLAQ 
surveys to form a survey of « 60000 QSOs. We find a hint of a peak in the QSO 
2-point correlation fun ction, £(s), at the same scale (« 105h _1 Mpc) as detected by 
lEisenstein et al.l f|2005T ) in their sample of DR5 LRGs but only at low statistical sig- 
nificance. We then compare these data with QSO mock catalogu es from the Hubble 
Volume N-body light-cone simulation used by Hovle et a l. (2002) and find that both 



routes give statistical error estimates that are consistent at « lOOh 1 Mpc scales. Mock 
catalogues are then used to estimate the nominal survey size needed for a 3-4er detec- 
tion of the BAO peak. We find that a redshift survey of « 250000 z < 2.2 QSOs is 
required over pa 3000deg 2 . This is further confirmed by static log- normal simulations 
where the BAO are clearly detectable in the QSO power spectrum and correlation 
function. The nominal survey would on its own produce the first detection of, for 
example, discontinuous dark energy evolution in the so far uncharted 1 < z < 2.2 
redshift range. We further find that a survey with « 50% higher QSO sky densities 
and 50% bigger area will give an w 6er BAO detection, leading to an error k, 60% of 
the size of the BOSS error on the dark energy evolution parameter, w a ■ 

Another important aim for a QSO survey is to place new limits on primordial non- 
Gaussianity at large scales. In particular, it is important to test tentative evidence we 
have found for the evolution of the linear form of the combined SDSS+2QZ+2SLAQ 
QSO £(s) at z w 1.6, which may be caused by the existence of non-Gaussian clustering 
features at high redshift. Such a QSO survey will also determine the gravitational 
growth rate at z ~ 1.6 via redshift-space distortions, allow lensing tomography via 
QSO magnification bias while also measuring the exact luminosity dependence of 
small-scale QSO clustering. 

Key words: quasar clustering 



1 INTRODUCTION 

Quasi-stellar objects (QSOs) have been used as tracers of 
large-scale structure for many years n ow. The first mea- 
surements were made by lOsmen |l981), then with the ar- 
rival of high -multiplex fibre systems, the sub j ect advanced 
rapidly (e.g. Boyle. Shanks. & Peterson 1988]; ICroom et al.l 
120051 ; Ida Angela et al.ll2008l ; iRoss et al.ll2009l and references 
therein) . Their clustering at small scales as measured by the 
correlation function is known to be consistent with the usual 
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7 = —1.8 power-law form for galaxies. The amplitude is 
comparable to galaxies at low redshifts and remains reason- 
ably constant with redshift. At larger scales the power spec- 
trum has been measured to be reasonably con sistent with 
the standard ACDM cos mological model (e.g. iHovle et al.l 
120021 ; lOutram et aT1l2005n . 

Here, we have combined the largest, spectroscopically 
confir med, QSO surveys from fibre spectrographs incl uding 
2QZ (ICroom et al.ll2004). 2SLAQ l|Croom et al.ll2009h and 
SDSS DR5 (jSchneider et al.ll2007l ) to form a redshift survey 
comprising some ?s 60000 QSOs. We have already found in 
these datasets that the small scale clustering of QSOs is 
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remarkably indepen dent of QSO luminosity at fixed redshift 
i|Shanks et alj|201ll ). 

In this paper, we first use the above combined survey 
to estimate the large-scale QSO correlation function and 
search for the BAO feature. We then outline our initial mo- 
tivation for an extended QSO redshift survey using the 'ef- 
fective survey volume' as a measure of clustering 'grasp'. 
We then make an empirical test of the errors on the QSO 
correlation function at « 100h~ 1 Mpc scales using both the 
data and mock QSO catalogues from the Hubble Volume 
N-body simulation. These routes allow a first estimate of 
the survey size needed for a significant BAO detection. We 
particularly focus on surveys that could be made wi th 2dF 
|Lewis et alj|2002h and AAOmega dSmith et alj|2004 ) at the 
AAT (c.f. the results of IWang et all (120091 1 for 3SO surveys 



s = 5.95 ± 0.45/i _1 Mpc and 7 
25h _1 Mpc range. 



1.16±2'i6 in the 1 < s < 



with LAMOST). We then use static log- normal simulations 
to test further BAO detectability and use Fisher matrix and 
Markov Chain Monte Carlo (MCMC) methods to test a 
QSO survey's competitiveness against other routes to the 
BAO scale and the evolution of w. Finally, we look at the 
prospects for detecting evidence for non-Gaussian cluster- 
ing at large scales via a QSO survey, in particular focussing 
on the possibility that the current QSO surveys show evi- 
dence for evolution in the linear regime of clustering, which 
might represent evidence for non-Gaussianity, if confirmed 
in a larger survey. 



2 QSO CLUSTERING DATA 

We start by making a study of the results from current QSO 
surveys (a) to see if the BAO peak can be detected in the 
correlation function and (b) to measure the errors to base 
an empirical estimate of the new survey size needed for an 
accurate BAO measurement. This survey size estimate will 
then be compared to those from simulations to determine the 
best QSO survey strategy in terms of area and magnitude 
limit. 



2.1 SDSS, 2QZ and 2SLAQ surveys 

Previously ICroom et all (|2005l ) used the 2QZ survey to es- 
timate the QSO correlation function and its dependence 
on redshift and luminosity. This survey contained ~ 22655 
QSOs in two w 375deg 2 NGC + SGC strips. The magnitude 
limit was 18.25 < bj < 20.85 and the re sulting QSO sky 
density was 31deg -2 . ICroom et al.l l|2005h measured so = 
5.4l^4g/i _1 Mpc and 7 = 1.2 ± 0.1 at 1 < s < 25h _1 Mpc for 
the amplitude and slope of the z-space correlation function, 
Ms). 

Ida Angela et al.l (|2008T ) then used the 2SLAQ survey 
of 9418 QSOs based on SDSS imaging to test the lumi- 
nosity dependence of the QSO clustering. The magnitude 
limit was 20.5 < qab{~ bj) < 21.85 and the resulting QSO 
sky density was ~ 45deg~ 2 , including the 2dF QSOs where 
the two surveys overla pped, in a total area of w 200deg 2 . 
Ida Angela et al.l (2008) found a £(s) amplitude and slope 
similar to that for 2QZ. 

Most recently, IRoss et all (|2009h have analysed the clus- 
tering of 30239 QSOs in the 4013deg 2 SDSS DR5 survey 
to iab = 19.1. This gives a QSO sky density o f 7.5deg -2 
in the uniform sample where Ross et all (|2009l ) measured 



2.2 Large-scale clustering comparison 

Here, the clustering analysis has been re-done to use consis- 
tent bins at comoving separations corresponding to BAO 
scales in the 2QZ, 2SLAQ and SDSS-DR5 spectroscopic 
QSO samples. The data and ran dom catalogues a r e the same 
as those used in t he analyses of Ida Angela et all (|2008T ) and 
iRoss et all (120091 1 for 2QZ+2SLAQ and SDSS-DR5, respec- 
ti vely. We used the 'UNIFORM' sample with 0.3 < z < 2.2 
of IRoss et all l|2009h which contains 30239 QSOs over 4013 
deg 2 . The 2QZ+2SLAQ sample contains 28790 0.3 < z < 
2.9 QSOs and its small- and inter mediate-scale clu s tering 
analyses have been pe r forme d by Ida Angela et all (|2008T l 
(see also ICroom et al.l (|2005l l for the 2QZ-only clustering 
analysis). We perform a new clustering analysis by count- 
ing pairs at separation, s, independently for the SDSS and 
2QZ+2SLAQ samples. The data-random, DR, and random- 
random, RR, pairs for each sample are normalised by N r d 
and N 2 d , respectively, where N r d is the ratio between num- 
bers of randoms and data (~ 20 for the 2QZ+2SLAQ and 
« 30 for the SDSS samples). The lLandv fc Szalavl |l993l l es- 
timator is then used to determine the £(s) from the summed 
pairs over the different samples. Note that the r e sults are in 
good agreement with those using the iHamiltonl (| 19931 ") esti- 
mator. 

In Fig. [T] we then compare the large-scale clustering 
results from the three surveys directly with each other us- 
ing the redshift-space correlation function, £(s). The cos- 
mology assumed in all cases is Q.a = 0.7, Q, m = 0.3. Fit- 
ting t he £(s) results consi stently in the 1 < s < 30h _1 Mpc 
range, Sh anks et al.l (|201ll l fitted real-space correlation func- 
tion scale-lengths, ro, assuming power-law slope, 7 = 1.8, 
infall parameter, /3 = 0.4, and line-of-sight pairwise veloc- 
ity dispersion < w 2 > J / 2 = 750kms _1 . These authors found 
agreement at the 1.4<r significance level in these results with 
SDSS giving r = 6.30±0.3h -1 Mpc, 2QZ giving r = 5.75± 
u.25h _1 Mpc and 2SLAQ giving r = 5.70 ± 0.35h _1 Mpc. 
The best overall fit is r = 5.90 ± 0.14h _1 Mpc. Thus the 
small-scale results suggest that it is reasonable to combine 
these 3 surveys and the £(s) result at large-scales is also 
shown in Fig. [1] The data is compared to a LCDM model 
and a 'wiggle-free' version of the model (|Eisenstein fc Hul 
1998). These models are normalised to match the data at 
s = 10 — 30h _1 Mpc (see Fig. [9] below). The errors are 
based on Poisson errors calibrated by jack-knife errors in 
all 3 cases. Generally the SDSS has the biggest errors, par- 
ticularly at the smaller scales. This is because of its rela- 
tively low sky density. This means that at « 100h _1 Mpc 
the 2QZ (and 2SLAQ) surveys dominate the statistics at 
the predicted BAO scale. The 2QZ survey error is ~ 2x 
smaller than the SDSS error and 1.8x smaller than the 
2SLAQ error. The blue filled circles then represent the over- 
all 2QZ+2SLAQ+SDSS result, produced by simply adding 
the QSO-QSO and QSO-random pairs across the surveys. 
The resulting error is fa 2 x larger than the error from the 
SDSS LRG sample of lEisenstein et all |2005l ). We see that 
there is some hint of a detection at 105h _1 Mpc but there 
is a similarly sized feature at m 85h _1 Mpc. Clearly while 
the data appear promising in terms of detecting the BAO 
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Figure 1. The large-scale redshift-space correlation function re- 
sults from the 2QZ, the 2QZ+2SLAQ and the combined SDSS 
DR5+2QZ+2SLAQ QSO surveys. The errors are empirically 
scaled by the average (1.2x) ratio of jack-knife to Poisson errors 
in this separation range. The results are compared to ACDM (red 
solid) and no-wiggle (green dashed) m odels and also the z = 0.35 
LRG result of lEisenstein et alj j2005l) . 

feature, a larger sample size is needed and in Section [5] be- 
low we will translate these empirical errors into a required 
survey size to measure the BAO scale. We shall also com- 
pare with the errors predicted by simulations and use these 
to optimise the properties of a new QSO redshift survey for 
large-scale clustering. 



3 QSO SURVEY EFFECTIVE VOLUME 

We next estimate the efficiency of future QSO clustering 
surveys via their effective volume (V e ff)- Although QSO sky 
densities are lower than for galaxy surveys, their clustering 
amplitude is relatively high, the exposure time required to 
establish redshifts is generally quite short, contamination 
rates are increasingly low and the volumes probed are very 
large. Moreover, the low QSO sky density can be viewed as 
an advantage in that it may well match the fibre density 
currently available from instruments like A AT 2dF. 

A rough measure of a survey's capability for power spec- 
trum or 2-point cor relation function analys is, the effective 
volume is defined by lEisenstein et alj {2005) and represents 
the survey volume that has a high enough QSO density for 
the shot noise to lie below the amplitude of a spatial power 
spectrum feature such as a BAO oscillation scale. The power 
spectrum or correlation function error at a particular scale 
then is proportional to V~fj 2 . In Fig. [2] we have calculated 
the effective volume for QSO surveys assuming the SDSS 
DR5 QSO n( z) in the redshift range 0.5 < z < 2.2 (see 
also Fig. 2 of I Wang et al . 2009). We have chosen a nomi- 
nal survey area of 3000deg 2 ; effective volumes of other sur- 
vey areas will scale linearly. Since the QSO n(z) is approx- 
imately independent of survey magnitude limit, the main 
other survey parameter is QSO sky density. We have cal- 



culated the effective volume at sky densities approximating 
those for the SDSS, 2QZ and 2SLAQ surveys at 10, 35 and 
80deg -2 . We also present the effective volumes at 140deg -2 , 
which is approximately the largest density accommodated 
by the 2dF fibre positioner (if tiling overlaps is considered, 
see later), and 280deg -2 which is the highest QSO density 
that is available from the Hubble Volume simulations (see 
Sect. 14.1] ). The assumed QSO correlation function amplitude 
was so = 6h~ Mpc which is also found to be alm ost indepen- 
dent of survey limit (see eg lShanks et al.ll201ll ). From Fig. 
[2] we see that QSO effective volume generally drops sharply 
as the spatial wavenumber increases; this drop is at a faster 
rate than for more highly sampled galaxy surveys such as 
WiggleZ. 

However, even at the 2SLAQ sky density of 80deg -2 , 
we see that at the scale of the first acoustic peak at k w 
0.02hMpc~ 1 , the effective volume of our nominal 3000deg 2 
survey overtakes that of the curre nt leading galaxy BAO 
survey, WiggleZ {Blake et al.ll201ll ). by a factor of « 3. Of 
course, even if the effective volume is only merely competi- 
tive with WiggleZ volume as it is at the 2nd and 3rd peak 
positions, then this still represents an advance, given the 
~ 3x higher redshifts of the QSOs than the WiggleZ galax- 
ies. To reach the same effective volume at the first acoustic 
peak of the current BOSS LRG survey at z « 0.55, a higher 
QSO sky density of 140deg -2 would be needed. This sky 
density would be reached at g sa 22. assuming a 10° 3m 
QSO number count slope {Boyle. Shanks, fc Petersonlll988l ). 
This QSO effective volume again applies at an ~ 3x higher 
redshift than the BOSS LRGs. Even the BigBOSS ELG sur- 
vey will produce an effective volume which is only ~ 2 x big- 
ger than for a 140deg -2 QSO survey when renormalised to 
the same area of sky. BigBOSS also has a significantly lower 
average redshift, z » 1. Thus the relatively crude effective 
volume measure suggests that a QSO survey of nominal area 
3000deg 2 and sky density in the range 80-140deg" 2 should 
produce large-scale clustering results that have similar pre- 
cision at the first acoustic peak to the BOSS LRG survey 
but at significantly higher redshift. 



4 BAO SEARCH IN SIMULATIONS 
4.1 Hubble Volume 

We next use the Hubble Volume simulation to measure the 
correlation function errors directly in mock surveys that 
range up to higher QSO sk y densities than 2QZ. This simu- 
lation {Evrard et al.ll2002i) used an initial mass power spec- 
trum with Q. b = 0.04, ficDM = 0.26, Qa = 0.7, H = 70 
km s _1 Mpc -1 and as = 0.9. Mo ck QSO cata l ogues from 
this simulation were generated by iHovle et all {20021 ). We 
did consider using newer simulations but although these fre- 
quently had higher resolution, they generally did not have 
sufficient volume to accommodate even the original 2QZ sur- 
vey. The Hubble Volume mocks are made in the form of past 
light cones as needed for accurate modelling of the QSO 
survey. The mock QSOs bias relativ e to the mass w a s mod - 
elled according to the algorithm of lHatton fc Colel {l998h . 

1 We assume that the QSO clusteri ng and bias continu e to be 
luminosity independent at this limit llShanks et alj|201l[) . 
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Figure 2. The effective volume as a function of wavenumber 
of spatial scale for a z ~ 1.5 QSO survey of area 3000deg 2 
and QSO sky densities varying from 10deg~ 2 to 280dcg~ 2 as 
shown. The assumed QSO correlation function amplitude was 
so = 6/i — ^^Mpc. Also shown is the effective volume for the BOSS 
LRG survey at z ~ 0.5 assuming sq = 10h -1 Mpc and lOOOOdeg 2 
area. The WiggleZ ELG survey at z ~ 0.5 has so = 4.4h _1 Mpc 
and lOOOdeg 2 area. The proposed BigBOSS ELG survey at z w 1 
has s = 4.4h _1 Mpc and 14000deg 2 area. The arrows indicate 
the positions of the 1st, 2nd, 3rd and 4th BAO peaks (from left 
to right). 

The main change from the previous 2QZ mocks is that the 
QSO sky density is approximately doubled from 35deg~ 2 
to 75deg -2 . The area of the mock survey is 15 x 75deg 2 or 
1.5 x the area of the 2QZ survey. Previously we have shown 
that our correlation function and power-spectrum estima- 
tion techniques can accurately retrieve the input functions 
in real and redshift space. The errors are jack-knife esti- 
mates based on 60 sub-samples from this contiguous area. 
The amplitude of the correlation function at small scales is 
ro = 6h _1 Mpc, close to the ro = 5.9h _1 Mpc shown by the 
data. 

Fig. [3] shows the large-scale correlation functions from 
the mocks with 11 25deg area and 80 and 280deg~ 2 QSO sky 
densities. We find that the errors are compatible with those 
extrapolated using simple Poisson scaling of the data-data 
pairs from the 2QZ+2SLAQ survey with its smaller QSO 
sky density and area. In fact, we find that in the relatively 
small area of the Hubble Volume simulation the BAO peak 
is barely detected, in either the mock at the standard QSO 
sky density of 80 deg~ or even at 280 deg" 2 . In the £(s) 
measured for the unbiased mass (not shown), the feature is 
detected but only at low significance, l-2(j. Thus although no 
feature is detectable in the relatively small Hubble Volume 
mocks, these data can still be used to estimate the likely 
errors in the w 3x bigger 3000deg 2 nominal QSO survey 
considered in Section [3] above. 

Fig. [4] shows the ratio of the jack-knife errrors (60 
sub-samples) from the above QSO mock catalogues from 
the Hubble Volume simulations. The ratio of the errors 
agrees with the Poisson prediction between the 35deg -2 and 
80deg -2 (also 105deg -2 see later) sky densities but the in- 
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Figure 3. Correlation function from QSO Hubble Volume mocks 
at 35, 80 and 280deg -2 compared to ACDM and no- wiggle models 
and also the observed data. The errors found for the 80deg -2 case 
have been scaled to a 3000deg 2 survey and applied to the ACDM 
model (red points + line). 

crease to 140deg -2 and 280deg -2 only achieves a factor of 
3.3 and 5 improvement in the error rather than the Poisson 
predicted factor of 4 and 8, respectively. So Poisson scaling 
works as far as the sky density of ~ 100deg~ 2 . 

4.2 Log-normal simulations 

We also r an static simulations s i milar t o the Gaussian simu- 
lation s of lBlake~fc Glazebrookl (|2003h ; I Glazebrook fe Blake 
(2005), drawing 3-D mode amplitudes according to power 
spectra for a standard ACDM model. The distribution used 
for these realisations was up-dated to log-normal rather than 
Gaussian to mimic better the effects of non-linearities in th e 
matter distribution (|Coles fc JoneJl99ll ; lBlake et alJl201.ll ). 
Fig- EJi shows the P(k) analysis of 400 simulations for a 3000 
deg 2 , z < 2.2 QSO survey with sky density 90deg~ 2 and a 
uniform space density of n — 1.6 x 10 -5 h 3 Mpc -3 . We split 
the redshift range into 3 parts, 0.4 < z < 1.0, 1.0 < z < 1.6 
and 1.6 < z < 2.2 with QSO sky densities of 18, 33 and 39 
deg -2 with an infall parameter of (3 = 0.58, 0.43, 0.32 and 
b = 1.4, 2.1, 3.0 in successive redshift ranges (|Croom et al.1 
2005). Although the sample is dominated by shot noise, 
BAOs are detectable in P(k) in the 2nd and 3rd slices with 
a precision comparable to WiggleZ and SDSS-LRGs. In Fig. 
[5^, the dotted lines are the result of an effective volume cal- 
culation for the errors, which agrees well with the scatter in 
the lognormal realizations. Th e accuracy in the P( k ) BAO 
comes from fitting the simple Ir31ake fc Glazebrookl |2003); 
I Blake et al.1 l|2006l ) model to the realizations. We detect the 
BAO in the 1 < z < 1.6 bin with ±5% accuracy for the 
BAO scale and in the 1.6 < z < 2.2 bin with ±3.7% accu- 
racy. Overall the BAO scale accuracy in the 1 < z < 2.2 
range and m 70deg -2 sky density is ±3%. 

Fig. [5p shows the mean correlation function result in- 
tegrated over the full 0.4 < z < 2.2 redshift range. Again 
we see the BAO feature clearly detected at w 3<r relative to 
zero signal. 
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5 QSO DARK ENERGY SURVEY 

5.1 Empirical survey parameter estimates 

We first take the empirical, jack-knife error estimates at 
100h -1 Mpc from the 2QZ £(s) result, also shown in Fig. Q] 
Ignoring the 2SLAQ survey simplifies the scaling of errors 
since in the Northern Cap at least, 2SLAQ relied on 2QZ to 
supply the QSOs with bj < 20.85. We can then simply scale 
the errors from 2QZ assuming a sky density of 35deg -2 and 
an area of 750deg 2 . The 2QZ amplitude and error in Fig. [T] 
suggest that a BAO peak at ~ 100ft _1 Mpc will appear at 
the m la level in a QSO sample of the current size (0.86a 
against the no- wiggle-model and 1.1a against zero correla- 
tion signal). Hence a « 4cr detection will require either a 16 x 
bigger survey at the 35deg -2 density or a 4x bigger survey 
at double the QSO sky density. The simulation results in Fig. 
[4] suggest that this Poisson sky density scaling continues at 
least as far as w lOOdeg -2 but not as far as 280deg~ 2 . Fig. 
also suggests that the error is 3x rather than 4x smaller at 
a QSO sky density of 140deg -2 compared to 35deg -2 , so a 
survey of only ~1300deg 2 would be required to achieve a 4a 
BAO detection. By mainly using previous 2QZ and 2SLAQ 
survey areas and taking tiling overlaps of « 20% into ac- 
count (to ensure survey areal completeness), the sky density 
for new QSO targets might only be « lOOdeg -2 , or « 300 
per 2dF field. If these QSOs could be efficiently detected 
with a contamination of only 25% or less then it may be pos- 
sible to achieve this density with only a single 2dF pointing 
per field. 2SLAQ achieved a 44% star contamination rate 
based on SDSS single epoch imaging data and improving on 
this rate mostly depends on achieving improved ugriz pho- 
tometry compared to SDSS. This should be possible using 
new surveys such as VST ATLAS (Shanks et al. 2011, in 
prep.). Note also that 2SLAQ used more traditional colou r 
cuts, and methods such as KDE (e.g. iRichards et al.l l2009). 
extreme-deconvolution etc., would improve target selection 
considerably. 
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Figure 4. The ratio of the jack-knife errors (60 sub-samples) 
from the QSO mock catalogues from the Hubble Volume N- 
Body simulation, using densities of 35deg — 2 , 80deg — 2 105deg — 2 , 
140deg" 2 and 280deg -2 . The final ratios of 2.22±0.28, 2.94±0.35, 
3.34±0.31 and 5.0±0.64 can be compared to the QSO sky density 
ratios of 2.29, 3.0, 4.0 and 8.0, confirming that the error scales as 
expected from Poisson statistics between 35 and as lOOdeg -2 . 



results of Fig. [T] again indicating that these survey param- 
eters will produce an ~ 4a detection of the BAO feature in 
f(s). Finally, the error on the £(s) peak at ~ 105h -1 Mpc 
in the log-normal simulations in Fig. [5p represents again an 
~ 3 — 4a detection. So there is generally excellent agreement 
between simulated estimates and the empirical results from 
Section[5T]that a 3000deg 2 QSO survey with an 80-90deg" 2 
QSO sky density will produce a significant detection of the 
BAO peak with a scale measurable to w ±3%. 



5.2 Simulated survey parameter estimates 



5.3 Comparison with other surveys 



We then take the jack-knife based error estimates from the 
Hubble Volume simulation. As noted above, we have checked 
that the jack-knife errors reduce approximately linearly as 
the QSO density increases up to « 80deg -2 . Fig. 2] shows 
that as the mock QSO sky density rose from 35 to 80, 140 
and then 280 deg -2 , factors of 2.3, 4 and 8.0, the jack-knife 
error in 60 subsamples reduced by factors of 2.22 ± 0.28, 
3.34 ± 0.31 and 5.0 ± 0.64. If we drop the z < 1 QSOs (see 
below) in the 140deg -1 mock catalogue, the sky density be- 
comes 105deg -1 and the error on the £(s) only increases by 
10%. Care was taken here that the small scale correlation 
function amplitude remained constant at the three sky den- 
sities so that the effect of sky density could be easily mea- 
sured. These ratios also reasonably agree with taking ratios 
of the average of the square root of the three lowest k effec- 
tive volumes in Fig. [5] So we again conclude that at least 
up to a sky density of ~ lOOdeg -2 , doubling the sampling 
rate approximately halves the error. We also see that for 
the 80deg -2 QSO mock, extrapolating the error from the 
mock survey area of 1125deg 2 to 3000deg 2 , indicates that 
the error is reduced to « 27% of the error in the 2QZ £(s) 



Fig. [6] compares the combined 1 < z < 2.2 BAO error 
for the dilation scale, Cy(z) = (Da(z)' 2 cz/ H(z)) 1 ^ 3 , of a 
nominal 3000deg 2 , 90deg -2 QSO survey with other cur- 
rent BAO surveys. The dilation scale is a measure which 
combines the information in the comoving angular diam- 
eter distance, Da(z), and the Hubble parameter, H(z). 
The errors a re no w generated using the fitting formula of 
iBlake et all (I2006T) which is calibra ted by lognormal real- 
izations (|Glazebrook fc Blake! 12005m for the BAO measure- 
ments. We see that the w ±4% error from the QSO survey 
is comparable to the BAO error at z = 0.35 from the SDSS 
LRG survey and WiggleZ at z = 0.6. Note that the error de- 
rived for the dilation scale from the QSO survey-specific log- 
normal simulations described above is ±3% which is clearly 
the more directly measured result. The results in Fig. [SI on 
the other hand, have the advantage that they are measured 
consistently between the various surveys. 

In response to a request from a referee, we note that 
the Hobby-Eb erly Telescope D ark Energy Experiment HET- 
DEX survey l|Hill et al l I200SI ) will allow fa 750000 Ly-a 
emitting galaxies to be mapped over 1.9 < z < 3.5 in 
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Simulations of a 3000 deg^ z<2.2 QSO survey 




Simulations of a 3000 deg^ z<2.2 QSO survey 



50 100 
s / h" 1 Mpc 
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Figure 5. (a) Predicted survey sensitivity for BAO from 400 log- 
normal simulations. Each successive set of data has been offset in 
y for clarity. QSO Power-spectrum BAO accuracy is 3% over full 
1 < z < 2.2 range in nominal survey, (b) The QSO £(s) from the 
log-normal simulations integrated over the 0.4 < z < 2.2 redshift 
range. The BAO signal is clearly detected. 



BAO forecast errors [Blake & Glazebrook 2006] 



o wiggleZ [measurement] 
• WiggleZ [forecast] 

o SDSS LRG [measurement] 
» SDSS LRG [forecast] 



I BOSS LRG [forecast] 



z<2.2 QSO [forecast] 



Figure 6. Predicted sensitivity for BAO in the nominal QSO 
survey consistently compared with other surveys. QSO power 
spectrum BAO a ccuracy forecasts from the fitting formula of 
iBlake et al] <2006l) . 



420deg 2 of sky in 150 clear nights. This survey is claimed 
to measure the BAO scale to « 1% accuracy. This survey 
will therefore produce similar errors to a QSO survey with 
an area of 4500deg 2 and a sky density of » lOOdeg -2 , in a 
mostly higher and hence complementary redshift range. 



5.4 Markov Chain Monte Carlo cosmology fits 

Finally, we ran some MCMC cosmology fits for WMAP dis- 
tance priors plus current and future BAO surveys for (f2 m , 
Q. m h 2 , wo, w a , f2fc) with the motivation of seeing if the data 
can distinguish curvature from evolving dark energy. For 
that reason we first focus on the joint likelihood of (w a , 
Slfe), which is shown in Fig. [7K for various combinations of 
surveys. We assume the CPL (jChevallier fc Polarskil lioOll ; 
iLinderlfeOQSt ) parameterisation for the evolving dark energy 
equation of state, i.e. w(a) — wo + w a (l — a). For the WMAP 
distance data we used the constraints on the shift parameter 
R, the aco ustic scale l a and the re dshift of recombination z* , 
as given in lKomatsu et all (|2009l ). 

Given just SDSS-LRG + WiggleZ (not shown), our 
nominal 3000deg 2 /80deg 2 QSO survey at z = 1.6 does 
give significant help. However, Fig. [7^. is based on WMAP 
+BOSS-LRG (+BOSS-Lya) and the accuracy of particu- 
larly the BOSS-LRG measurement at z ~ 0.5 provided by 
this combination means that the QSOs (green contour) will 
only decrease the errors in w a , fife by at most 10-20% in this 
parameterization. Clearly in the case where there is little 
evolution and w a « 0, even in their more restricted redshift 
ranges, BOSS LRG and Ly-a surveys already constrain w a 
as strongly as the nominal 3000deg 2 +80deg -2 QSO survey 
at z ~ 1.6. The same result holds in the Wq, w a plane in Fig. 
[7p. On the other hand, it should be noted that the nominal 
QSO survey would on its own produce the first detection of, 
for example, discontinuous dark energy evolution in the so 
far unexplored 1 < z < 2.2 redshift range. 

We next consider what survey parameters would lead to 
significant improvements in the cosmological forecasts, even 
in the case where the dark energy evolves smoothly from 
z ~ 0.5 to z » 1.6. Bigger QSO surveys shown in Fig. [7^,b 
assume 1% and 2% errors in Dv and the contours suggest 
that a 1.5% error is likely to give significantly smaller errors 
than the competing BOSS LRG+Lyman-a surveys in the 
w a ,ftk and the wo,w a planes. 

If so, then we would need to increase the QSO sky den- 
sity by a factor of ~ 1.6 to ~ 140deg -2 , reducing the QSO 
Dy error by a factor of ~ 1.5 from 3% to 2%. This would 
require a 0.6mag fainter mag limit taking us to g < 22.5 
rather than g < 21.85. Then we could drop the z < 1 QSOs 
which takes us back to ~ 105 deg -2 (see Section 5.2), with- 
out much loss of BAO S/N as evidenced from the log-normal 
simulations and Fig. 4. This should be possible with VST 
ATLAS imaging data. If we wanted to get below a 2% error 
then an additional « 50% of the area ie « 4500deg 2 would 
give a BAO error of 1.6%. This may be possible in a 1- 
2hr AAOmega exposure time and the survey would then be 
completed in ~ 200 clear nights. From Figs. 7, such a QSO 
survey would produce an error on the dark energy evolution 
parameter, w a , which is ~ 60% the size of that from BOSS 
LRGs. The BAO detection significance would be ~ 6a as 
opposed to 3 — 4a in the nominal 3000deg 2 , 80deg -2 QSO 
survey. 
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Figure 7. Predicted survey constraints for {1%, m a , mio from BAO 
MCMC fits, lcr contours and errors are shown, (a) Our nomi- 
nal QSO survey produces a 3% BAO error (green contour)and 
will provide the first detection of any non-standard (eg discon- 
tinuous) evolution in the dark energy equation of state in the 
1 < z < 2 range. But if the evolution remains relatively stan- 
dard/continuous then a survey with a 50% bigger area and a 50% 
higher sky density will produce a 1-2% QSO BAO error which 
will significantly improve over BOSS constraints on w a and 
(cyan, blue contours), (b) Similarly, even in the case of stan- 
dard/continuous dark energy evolution, a survey with a 1-2% 
BAO error will provide significantly improved constraints over 
BOSS in the (wo,w a ) plane (cyan, blue contours). 



5.5 Redshift-space distortions 

The nominal 3000deg 2 redshift survey of 250000 QSOs 
would also be able to probe cosmology via redshift-space dis- 
tortions. Now a redshift-space distortion test of non-Einstein 
gravity is more difficult at high redshift because Q m (z), at 
least in FLRW-based models, tends to unity at high z , mak- 
ing t he 7 index of the gravitational growth rate |Linderl 
120051 ). f{z) = f2 m (z) 7 , more difficult to determine. However, 
interesting cosmological constraints can still be obtained. 
The infall parameter governing redshift-space distortions, is 
defined as j3 = Q.^ n /b, and so depends on the bias, 6, as well 
as the gravitational growth rate, fiJJ,. Previously for the 2QZ 
and 2SLAQ QSO surveys we have used redshift-space dis- 
tortion and the evolution of the QSO clustering amplitude 
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Fishei\ matrix forecast growth errors [White et al. 2009] 

\ 

^ o WiggleZ [measurement] 




z<2.2 QSO [forecast] 



BOSS tRG [forecast] 



Figure 8. Predicted sensitivity for growth of structure in the 
nominal QSO survey consistently compared with other surveys. 
Growth of stru cture forecasts are genera t ed us ing the Fisher ma- 
trix formula of lWhite, Song, fc Percivall ll2009h . 



to solve for Q m (z = 0) and bias, b(z = 1.6) simultaneously 
(da Angel a et al.l2005l ; lda Angela. Outram. fc Shanks 2005; 
Ida Angela et al.ll200&t ). This test also involves the Alcock- 
Pacynzski geometric test. The bias, b(z = 1.6) can then be 
used to derive the amplitude of mass clustering, ie as, at 
2 = 1.4. 

Recently, the combination / x <r™ ( = /3 x <j| if b = 
a^/a™) has become the prime target for redshift-space dis- 
tortion studies, since it can discriminate between modi- 
fied gravity models wi thout needing to determine the bias 
l|Song fc Percivalll2009l ). Redshift space distortions can thus 
provide a strong test of Einstein's gravity independently 
of geometrical cosmologi cal tests that use s tandard candles 
and rods, such as BAO (|Guzzo et al.ll200^ ). Redshift-space 
distortions can further be used to give an estimate of the 
masses and hence ma ss-to- light ratios of galaxy group haloes 
in CDM models (e.g. iMountrichas et alll2009l ). 

In Fig. [HI we show the error in the gravitational growth 
rate - mass fluctuation parameter, fas, for the QSO survey 
as estimated using the publ ic ly ava ilable Fisher matrix code 
of IWhite. Song, fc Percrval |2009l ). The result is averaged 
over the full range of clustering scales 1 < s < 100h _1 Mpc. 
We see th at the overall resu lt is again comparable to that for 
WiggleZ (|Blake et a,l.ll201 lh and the SDSS LRG survey but 
at a significantly higher redshift. If we made a survey at the 
~ 140deg -2 QSO density and over an area of 4500deg 2 , the 
error on these redshift space distortion results would reduce 
by a further factor of ~ 1.8x. So with these parameters the 
error on fas would reduce from 6% to 3.3% in this redshift 
range making the survey even more competitive in the 1 < 
z < 2.2 range. 



6 QSO CLUSTERING AS A PROBE OF 
NON-GAUSSIANITY 

6.1 Background 

Inflationary models with standard slow-roll inflation will 
produce very little non-Gaussianity, but models that de- 
viate from the slow-roll assumption can produce a sig- 
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nificantly non-Gaussian seed field (see iBartolo et alj |2004| 
for a review). Non-Gaussianity is normally parameterised 
through some amplitude of a quadratic term in the pri- 
mordial Bardeen potential, Jnl- This parameter gives the 
coupling between a triangle of three k-modes, and can ei- 
ther be applied to the local ('squeezed' isosceles triangles) 
or equilateral form of non-Gaussianity. It can be measured 
using the large scale galaxy Bispectrum. The best current 
limits on non-Gaussianity are obtaine d from the cosmic mi - 
crowave background radiation (CMB; iKomatsu et al.|[201ll ) 
who found /^£ ai = 32 ± 21. Galaxy surveys can actu- 
ally be more powerful as they can sample more modes in 
3-D space than the CMB can on the 2-D surface of the 
sphere. It also samples the structure of matter perturba- 
tions at smaller scales, making it complementary to ex- 
periments such as Planck. By probing scales in the range 
k = 0.01 — 0.2/iMpc - , it will be able to l ink constraints 
from the CMB with constraints from clusters l|Lo Verde et al.l 
2008), constraining possible scale-dependence of the non- 
Gaussianity. 

In galaxy surveys, non-Gaussianity can produce a 
scale-dependent boost of the halo power-spectrum at k < 
O.OSftMpc^ 1 and this evolves as (1 + z). Although this can 
potentially be confused with the full general relativistic cor- 
rection of the galaxy power spectrum at k $J O. OlfaMpc , 
the effect becomes important only beyond z w 3 (|Yodl2010l) . 
Hence QSO surveys with their large volumes and redshift 
range s make an i deal ba sis for this test (e.g. ISlosar et al.l 
120081 ). IXia et all (|2010al lbl) has recently produced upper 
limits on non- Gaussianity from th e SDSS-DR6 photo-z 
QSO catalogue jRichards et al.ll2009l) an d The NRAO VLA 
Sky Survey (NVSS; ICondon et all 1 1993 ). The NVSS auto- 
correlation function sho ws some evidence for a positive tail 
extending to 5-6 degrees et al. 2010a, confirming previ- 
ous results from lBlake fc Wallll2002r ) which could be caused 
by non -Gaussianity, implying f^l al = 62 ± 27. IXia et al.1 
|2010bl ) found lower but still consistent angular correlation 
functions from a million QSOs in the SDSS DR6 dataset, 
implying /^£ ai = 58 ± 24. However, /jvx measurements 
from high-z photometric surveys can contai n systematic bias 
due t o gravitational lensing magnification (jNamikawa et all 
2011). QSO redshift surveys will provide more accurate and 
stringent constraints. 

6.2 Non-Gaussianity constraints via the nominal 
QSO survey 

A 3000deg 2 QSO survey will give highly competitive con- 
straints on primordia l non- Gaussianity in the density field. 
ISefusatti fc Komatsul |2007l ) calculated how effective future 
galaxy surveys would be at measuring the /ml parameters, 
simultaneously with the non-linear bias. Their Fig. 6 shows 
predictions of the senstivity of surveys of different volumes 
with a galaxy density of 5 x 10~ 4 (h/Mpc) 3 . Their fore- 
casts demonstrate that our nominal QSO survey has cer- 
tain advantages over its low redshift counterparts. By mak- 
ing measurements at higher redshift, it is less affected by 
non-linear structure formation, and can measure the Bis- 
pectrum down to smaller scales. Its main advantage will 
be its volume, which will be larger than most other funded 
surveys planned for z > 1. We estimate that our survey 
will be able to constrain /jvl (local) with an error of about 
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Figure 9. The combined DR5+2QZ+2SLAQ QSO correlation 
function from Figs. [T]now plotted on log axes and compared to 
the 2dFGRS £(s) at 2 = 0.12 of Hawkins et al and the scaled 
SDSS LRG of Eisenstein et al and a linear model fitted to 
LRG surveys in the range 0.35 < z < 0.7. The slope of the QSO 
£(s) at z fa 1.6 is marginally (2<r) flatter than the linear model 
fitted at z < 0.7. 



±15, and /nl (equila t eral) o f about ±150, according to the 
ISefusatti fc Komatsul (|2007l ) analysis. This is better than 
any current survey : these authors predict uncertainties al- 
most 10 x larger for the SDSS-LRG survey, for example. Our 
estimate is similar to the best current CMB result, but this 
would be the first competitive test made using QSOs. 



6.3 Testing for non-Gaussianity via large-scale 
clustering evolution 

Recently, angular c orrelation functio n studies of EROs in 
SA22 at z « 1.5 jKim et al.1 l201ll) and LRGs in SDSS 
Stripe 82 at z fa 1 (fNikolouda kis et al.ll201l] ) have shown ev- 
idence for a flatter slope in the range 10 < r < 100h _1 Mpc 
compared to the lower redshift, SDSS {z fa 0.35), 2SLAQ 
(z fa 0.55) and AAOmega LRG (z w 0.68) surveys. The evo- 
lution is small but statistically reasonably significant (fa 3er). 
Since evolution is not expected in the linear regime in the 
standard cosmological model, one interpretation of this evo- 
lution is that it might correspond to evidence for a non- 
Gaussian feature in the galaxy correlation function becom- 
ing more prominent at hig h redshifts, simi l ar to that found 
above bv lXia et all (|2010ai ); [Blake fc Walll (|2002l ). 

However, the angular correlation function is particularly 
susceptible to artefacts in the data such as small artifical gra- 
dients. The amplitude of the redshift space correlation func- 
tion is intrinsically higher because of the lack of projection 
effects. Therefore it is interesting to use the QSO clustering 
correlation function, £(s), to look for evolution in the large- 
scale slope of the correlation function. In Fig. [9] we therefore 
compare the combined QSO correlation func tion to the 2dF- 
GRS f(s) at z « 0.12 (lHawkins et al.l|2003l) and the SDSS 



2005). We also show 



LRG £(s) at z fa 0.35 l|Eisenstein et al 
a linear model that was fitted to correlation functions from 
1.5 million LRGs in three pho tometric samples with a verage 
redshifts z = 0.35,0.55,0.68 (jSawangwit et alj|20lj ). This 
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model assumes a CDM Universe with = 0.73, fi m = 0.27, 
fbaryon = 0.167, a 8 = 0.8, h = 0.7 and n s = 0.95. This 
model has also been corr ected for scale-dependen t redshift- 
space distortion following lEisenstein et all |2005l ). It can be 
seen that this model also gives an excellent fit to the 2dF- 
GRS [z = 0.12) and SDSS LRG (z = 0.35) 3-D correlation 
functions, £(s), as shown. At the level of the current errors in 
the observed correlation functions, it is not possible statis- 
tically to distinguish the form of the observed low-2 galaxy 
and high-2 QSO correlation functions. However, when \ 2 
fitted in the 15 < s < 95h~ 1 Mpc range, the slope of the 
z ~ 1.6 QSO correlation function appears flatter at the 
2a level than the linear model fitted to the lower redshift 
(z < 0.7) surveys^ This provides some limited support to 
the results from the LRG angular correlation functions at 
8 = 1 and z — 1.5 but clearly more data are needed. Even 
the nominal QSO survey would provide an ~ 3 — 4x reduc- 
tion in the errors in the 10-100h _1 Mpc range and allow a 
much more significant search for non-Gaussian evolution at 
large QSO separations. Furthermore, the possible detection 
of non-Gaussian evolution in the current LRG and QSO sur- 
veys is a potent reminder that dark energy evolution at high 
redshift may also be of an unexpected form which should be 
measured rather than assumed. 



7 OTHER QSO Z SURVEY SCIENCE 

More QSOs will also provide new data on the small-scale 
amplitudes of QSO clustering as a function of luminosity 
and will improve £(s) statistics at SDSS, 2QZ and 2SLAQ 
depths and this has importance for QSO formation and 
evolution models. The amplitude of QSO clu stering seems 



rema rkably independent of QSO luminosity dShanks et al 
12011]). is marginally significant (|Porciani fe Norberd 1200*6 



and Shcn et al 



2009), and so checking for any small luminos- 



ity dependence can put strict limits on models of QSO envi- 
ronment and their host halo and BH mass. A 3000deg 2 sur- 
vey will contain 30000 QSOs at the SDSS depth, « 100000 
at the 2QZ flux limit and « 120000 at the 2SLAQ limit, 
so there will be significant improvements in the small-scale 
QSO clustering measurements in all 3 luminosity ranges. 

A further application of a QSO survey is QSO 
lensing via magnification bias llMvers et al. I 120031 . 120051 ; 
IScranton et al.ll2005l ; lMenard et alj|20ld ). Myers et al found 
stronger results using QSO s pectroscopic z than Me n ard et 
al did using QSO photo-z. iMountrichas fc Shanks] (|2007m 
found that contamination of the photo-z QSO sample by 
low redshift QSOs could reduce the anti-correlation signal 
that is expected from lensing of faint QSOs. A large spec- 
troscopic QSO survey would reduce the errors on the 2QZ 
results significantly and test the validity of the photo-z re- 
sults, as well as complementing galaxy shear weak-lensing 
tomography. 

QSO surveys also allow investigations of the topol- 
ogy of the Unive r se. Fo r example, using the 2QZ survey, 
IWeatherlev et all (|2003h searched for exact QSO spectral 
pairs at large separations to check for topologically closed 



2 Since the jack-knife errors are approximately Poisson in this 
range, we have ignored the covariance between points in our x 2 



analysis. 



universes. Again a larger QSO survey would allow more 
stringent constraints on such models. 



8 SUMMARY AND CONCLUSIONS 

We first made a new correlation function analysis of the 
combined SDSS, 2QZ and 2SLAQ QSO surveys, compris- 
ing some 60000 QSOs. We focussed on the large-scale, 
s ~ 100h _1 Mpc, results to test the strength of BAO sig- 
nal that could be detected in the current dataset. We found 
that the «22000 2QZ QSOs dominate the signal; although 
2SLAQ has a higher sky density it has only «6000 QSOs 
in total and although SDSS has a larger number of QSOs 
its contribution is less significant because of its low sky den- 
sity. We observe a possible peak at » 105h~ 1 Mpc wh e re the 
SDSS LRG £(s) peak was found by lEisenstein et all (|2005h 
but here it is only detected at a low significance of ~ lo- 
in the combined dataset and other peaks are seen at other 
separations at a similar significance. 

We then proceeded to investigate the QSO survey pa- 
rameters that would be needed to make an « 4<r detection 
of the BAO peak. We conclude that our nominal survey of 
250000 QSOs in a sky area of 3000deg 2 will allow us to make 
a 4<r detection of the BAO scale at 1 < z < 2.2, an as yet 
unexplored range for cosmology. This ±3% BAO scale mea- 
surement will determine the high-redshift evolution of the 
dark energy equation of state, p — w(z)p, and in particular 
show if there is any large (> 15%) deviation from w — — 1 
in the 1 < z < 2.2 redshift range. But even if w a ~ 0, then a 
survey with 50% higher QSO sky density and a 50% bigger 
area will approximately halve the BAO error to 1.6%. At 
this point the QSO survey will also approximately halve the 
error on the dark energy evolution parameter, w a , and sig- 
nificantly reduce the errors on wo and fij, compared to the 
BOSS LRG and Lyman-c* BAO results at lower and higher 
redshifts. 

A QSO survey can also set powerful new limits on the 
existence of non-Gaussianity (/nl) in the primordial density 
field. We have found that the combined SDSS+2QZ+2SLAQ 
QSO survey shows possible « 2a evidence for evolution in 
the linear form of the z ~ 1.6 QSO correlation function, in 
the sense that it shows a flatter slope than a linear model 
fitted to galaxy surveys at z < 0.7. Even the nominal QSO 
survey would improve the significance of detection of this 
slope change by 3 — 4x. If confirmed, the result could in- 
dicate the detection of the evolution of a non-Gaussian fea- 
ture in the large-scale QSO correlation function. The tenta- 
tive evidence found here and in projected galaxy correlation 
functions for unexpected evolution of large-scale structure 
emphasises that dark energy evolution at high redshift may 
also be of an unexpected form. 

The QSO survey will further support analyses of red- 
shift space distortions to measure f x as and test Ein- 
stein gravity versus other modified gravity models. We 
shall also make the mo st rigorous application of the 
lAlcock fc Paczvnskil (|l979l ') test so far of the prediction of 
the standard ACDM model for the amplitude of the mass 
clustering, as, at z » 1.5. Other science it will do includes 
making the most accurate determination of the luminosity 
dependence of QSO clustering at small scales in order to 
probe QSO formation and evolution via QSO environment. 
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The survey could also use QSO lensing magnification bias 
to measure the mass and the bias of foreground groups and 
clusters and to do lensing tomography. 
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